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OF INDIVIDUAL RANDOM NUMBERS OF A 
SEQUENCE OF RANDOM NUMBERS OF A 
1/F NOISE 

5 

CROSS-REFERENCE TO RELATED 
APPLICATIONS 

This application is a continuation of copending Interna- 
tional Application No. PCT/DEO 1/04376, filed Nov. 22, 10 
2001, which designated the United States and was not pub- 
lished in English. 



BACKGROUND OF THE INVENTION 
Field of the Invention 



15 



The invention relates to a method for generating sequences 
of random numbers of a 1/f noise. 

Random numbers of a 1/f noise can be used, for example, 20 
in a transient circuit simulation that takes account of noise 
influences. The 1 /f noise is understood as a stochastic process 
with a specific frequency spectrum that can be described by 
the equation 

25 



of the random numbers are determined, and this is used for a 
simple and accurate generation of the corresponding random 
number sequences. 

The method according to the invention for generating at 
least one sequence of random numbers of 1/f noise first pro- 
vides the following steps in this case: 

a) determining a desired spectral value P; and 

b) determining an intensity constant const. 
Characteristics of the 1/f noise to be simulated are thereby 

axed. 

As a result, the number of random numbers of the 1/f noise 
that are to be generated, and a starting value for a running 
variable n used for the simulation are fixed. 

Until the desired number of elements y(n) of one or more 
vectors y of length n are calculated from 1/f-distributed ran- 
dom numbers, the invention provides the loop- type repetition 
of the following steps: 

a) increasing a current value of a running variable n by 1; 

b) fixing a simulation time step [t„_j ; t J ; 

c) determining an element C, y of a covariance matrix C of 
dimension (nxn) according to: 



:=const.(-\tj- ti f +l +1/;-, -at +1 + \t } - - -h-xt* 1 * 



i, j = 1 n 



30 



1/f noise sources are suitable for modeling noise influences in 
a multiplicity of technical and physical systems, and for sys- 
tems for estimating and forecasting events on the financial 
markets. In particular, many electronic components such as, 
for example, pn diodes and MOS field-effect transistors 35 
exhibit 1/f noise sources. 

It is possible to approximate 1/f noise sources by carrying 
out a summation of the effects of many noise sources which 
^each exhibit a Lorentz- spectrum as a frequency spectrum. 
Such noise sources can be modeled, for example, by the 40 
system response of a linear time-invariant system, which can 
also be denoted as an LTI-system, at whose system input 
white noise is applied. It is disadvantageous in this mode of 
procedure for the dimension of the system of differential 
equations that are solved numerically to be excessively 45 
inflated. The result of this is long computing times and a high 
storage requirement of a computer system that is used to 
simulate a system that is subject to the influence of 1/f noise. 



SUMMARY OF THE INVENTION 



d) determining a matrix C" 1 by inverting the covariance 
matrix C; 

e) determining a variable a in accordance with the rule 

G=sqrt(\/e(n,n)\ 

where sqrt denotes the "square root" function, and e(n,n) 
denoting the element of the inverted covariance matrix C 1 
indexed by (n,n); 

f) determming a (0,l)-normally distributed random num- 
ber which forms the nth component of a vector x of 
length n; 

g) forming a variable \i from the first (n-1) components of 
the nth row of the inverted covariance matrixC" 1 and the 
(n-1) elements of the vector y. which were calculated for 
a preceding (n-1) simulation time step, specifically in 
accordance with the following rule: 



It is accordingly an object of the invention to provide a 
method for on-demand generation of individual random num- 
bers of a sequence of random numbers of 1/f noise that over- 
comes the above-mentioned disadvantages of the prior art 55 
methods of this general type, which can be carried out quickly 
and with a low computational outlay. It is also an object of the 
invention to specify an improved method for simulating a 
technical system that is subject to 1/f noise. Finally, the aim is 
also to specify a computer system having a computer program so 
for determining sequences of random numbers of 1/f noise 
that can be executed quickly and that lays claim to only a few 
resources of a computer system. 

In accordance with the invention, the problem of noise 
simulation in the modeling of the system to be simulated is 65 
transformed into the problem of generating a random number 
sequence. In accordance with the invention, the correlations 



50 



y ( „_ 0 denoting the first (n-1) components of the vector 

y» 

C^ 1 denoting the first (n-1) components of the nth row 

a of the inverted covariance matrix C" 1 , and C^, denot- 
ing the element of the inverted covariance matrix C 1 
indexed by (ri,n); and 
h) calculating an element y(n) of the vector y of length n 

from 1/f-distributed random numbers according to the 

following rule: 

Simulations of technical systems can be arbitrarily length- 
ened with the aid of the method according to the invention. It 
is possible for this purpose to generate additional 1/f-distrib- 
uted random numbers in a simple way when already gener- 
ated 1/f-distributed random numbers are to hand. Moreover, a 
simulation can be based on results of previously simulated 
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time intervals. This so-called restart capability constitutes a e) forming q variables u* in accordance with the -following 
very important property for the practice of simulation. It is rule: 
precisely for 1/f noise sources that this can be achieved only 
with difficulty, because random numbers that simulate a 1/f 

noise source for a specific time interval depend on random 5 y£-nr£^ 
numbers for earlier time intervals that have already been w : = 

determined numerically. The present invention also permits 
the use of an adaptive step size control without thereby sig- 
nificantly increasing the computation times for the simulation 

of the technical system. Such an adaptive step size control 10 y (n _ nJc denotes the first (n-1) components ofthe vector^ that 
considerably increases the precision and the computation were calculated for a preceding simulation time step. C. t „ 
time efficiency in the numerical determination of the dynamic denotes the first (n-1) components ofthe nth row ofthe 
performance of a simulated technical system. inverted covariance matrix C" 1 and denotes the element 

It is no longer necessary in the case ofthe method accord- Qf ^ inverted covariance matrix CT 1 indexed by (n,n). This is 

ing to the invention to prescribe the time interval to be simu- 15 • d t f()r k=1 q 
lated. Precisely by the provision of variable step sizes, it is 9 ' ' * ' 4 " 

also possible to perform an adaptation to current system f) calculating q elements y M which form the respectiventii 
dynamics, and this increases the accuracy ofthe simulations. component ofthe vector y* of length n from 1/i-distab- 

The present invention specifies a method for generating uted random numbers, specifically according to the tol- 

sequences of 1/f-distributed random numbers successively, 20 lowing rule: 
that is to say element by element. In this case, the method ^ # 

ensures that each newly generated random number depends in 
a correct way in the stochastic sense on the previously gen- 
erated 1 /f-distributed random numbers . It is thereby possible k-i , . . . , q. 

to generate me respectively required random numbers in the 25 ^ vectors (k=i } . . . , q) of length n from 1/f- 
course ofthe numerical simulation of a circuit. distributed random numbers are arranged with particular 

The invention uses the theory of conditional probability advantage m a NOISE matrix and specify in a simulation the 
density in order to generate a 1/f-distnbuted random number influences of a system to be simulated, 

that correctly ensures the stochastic relationship ofthe ran- . , - _ • lllat 

dom number with the already generated random numbers 30 In accordance with the invention, the concept for Simula - 
required for preceding simulation steps. ing 1/f noise is based according to the invention on the tol- 

In a particularly advantageous refinement of the method lowing train of thought. The dynamic performance of a sys- 
according to the invention, q sequences of random numbers of tem mat is exposed to stochastic influences is modeled 
a 1/f noise are calculated simultaneously, instead ofthe steps adequately by a stochastic process. In order to simulate such 

to be repeated in the loop-type fashion: 35 a system dynamics, the general approach is to calculate 

length n; of ^ fundamental stochastic processes. Sunulating systems 

b) forming'a variable u from the first (n-1) components of with 1/f noise sources requires paths of stochastic integrals of 

menmrowofmemvertedcovariancematrixC'^ndthe 40 the form 
(n- 1 ) elements ofthe vector y_ which were calculated for 



a preceding (n-1) simulation time step, specifically in 
accordance with the following rule: 



45 



Jo 7 



(s)ds 



f i: = 7^1 



(s)ds to be calculated numerically. Here, s is an integration 
variable and t (upper integration limit) denotes the time, 

50 

where y (/1 _l) denotes the first (n-1) components ofthe vector 

y_, Cr* denotes me first (n-1) components ofthe nth row of m{s)ds 

meTinverted covariance matrix C"\ and C^, denotes the 7 

element of the inverted covariance matrix CT 1 indexed by 

(n,n); and 55 

c) calculating an element y(n) of a vector y of length n from ( s )ds denotes a 1/f noise source, andY(s) a stochastic process 
1/f-distributed random numbers according to the follow- that describes the temporal dynamics of a variable, for 
ing rule: example the electric voltage in the circuit simulation. 

Y(n)=x(n)*<HK Using B FB Js) to denote that stochastic process whose 

60 derivative (mathematically: derivative in the distribution 
the following steps are provided: sense) yields the 1/f noise process 

d) determining a quantity q of (0,l)-normally distributed 
random numbers x M which form the respective last 
component ofthe vectors x* of length n, k=l , . . . , q. It is J]i (5) 
to be borne in mind in this case that the respective first 65 7 
(n-1) components ofthe vectors x* have already been 

calculated previously in the step. 
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it is possible to write the stochastic integral to be calculated as 



f t Y(s)rj l (s)ds= P Y(s)dB FBM (s) 
Jo 7 Jo 



(1.1) 



The integral on the right-hand side is to be interpreted as a 
Riemann-Stieltjes integral of the stochastic process Y(s) with 
the process B^^s) as an integrator. This integral can be 10 
approximated by a sum by decomposing the integration inter- 
val [0,t] in accordance with 0sto<tj< . . . <m into n disjoint 
subintervals [t,., t,_ y ], i=l, n: 



Since the individual random variables &B FBj Ji) represent 
increments of a Gaussian stochastic process, the random vari- 
able vector (AB FU J}1 • • • AW n » determines an n-di- 
mensional Gaussian-distributed random variable, and thus 
completely defined by its (n-dimensional) expectation E and 
its covariance matrix C. The two variables can be calculated 
from the formulas (1 .3) and (1 .4) as 



C :=Cov(AB f B«(i).^ fSM (j)) = 
— w 



(3.5) 



(3.6) 



rt n 
Jo [illegible] 



(1.2) 
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This sum is a random variable. The dependence on the 20 
result co of the random experiment has been consistently 
omitted 

A process B Fm £s) whose generalized derivative has a 1/f 
spectrum is known in the literature under the name of 'Frac- 
tional Brownian Motion'. B Fm Jj>) is a Gaussian stochastic 25 
process and is completely characterized as such by its expec- 
tation 



and by its covariance function 
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The method according to the invention for on-demand 
generation of suitable random numbers essentially reduces 
the simulation of 1/f noise influences to the generation of 
realizations of the random variables PraJt/HiWt-i)]* 
that is to say to the generation of increments of Fractional 
Brownian Motion. 

The present invention permits the online generation of the 4Q 
required realizations of the random variables: &B FBA Ai), that 
is to say in the course of the successive integration of the 
system equations. This results in two requirements placed on 
the method: 

(a) The length n of the sequence of random numbers 45 
{ABraXl), . . . ,AB/rgaX n )} m^t remain variable dur- 
ing a simulation run. In particular, it must be possible at 
any time to lengthen the simulation (restart capability). 
This implies the capability of the method to generate the 

. .additional random numbers required. for this purpose 50 
such that they are correlated in a correct way with the 
subsequence already generated. 

(b) Let ti be the currently attained time in the course of a 
simulation. It must then be possible to determine the 
time interval [t,., t /+1 ], that is to say the step size of the 55 
next integration step, from the instantaneous system 
dynamics — that is to say adaptively. 

The invention meets both requirements by specifying a rule 
as to how a realization of {AB^^l), . . • ,AB F5A .Xn)}, that is 
to say a sequence of random numbers, can be generated 60 
successively, that is to say element by element. The step size 
At,:=t -t,_i can be freely selected in this case for each new 
random number. 

The first step is to investigate the approach for so-called 
"conditional densities". 65 

The first step is to consider the distribution of the random 
variable vectors (AB^^l), . . . .AB^^n)). 



const- (-If,- - Ut 1 + - «f* 1 + 

|/;W^ +1 H^i-'Ml^W = l n 

The online method according to the invention is now to be 
specified in the form of a complete induction. 

The basis of the induction— and thus the starting point of 
the method— is the realization of a real Gaussian distribution 
with expectation 0 and variance Z n =2 const ! At, l^ 1 . 

It must be specified for the purpose of inductive influence 
how to expand a realization of (AB F57W 0)> • • • .AB/W 11 - 1 )) 
by a realization of AB^/n) so as to produce overall a real- 
ization of (ABjWD, • . . ,AB Ff?JW (n)). In order to simplify the 
notation, the already "cubed" subsequence of random num- 
bers may be denoted by (y u . . . y„_i)=:y<„_n r and the real- 
ization, still to be cubed, of AB^^n) may be denoted by y„. 

The problem can now be formulated as now described. 

Let there be given an n-dimensional zero mean Gaussian 
random variable Z with the covariance matrix C . Let me first 
n-1 elements of a realization of Z be already cubed and 
known in the form of a random number vector y^iy 

The aim is now to find the distribution from which the nth 
element y n must be taken in order to complete } to form 
a realization y.=(y (n _i) > yj of Z. 

A solution to this task can be found when the conditional 
probability density ;f(y n ly (fI _i)) are considered for y— with 
the condition that y (n _ n already holds. This variable can be 
calculated in the present case of a Gaussian normal distribu- 
tion as 



In 



(3.7) 



The variable C^ n is yielded in this case from the following 
notation of the inverted covariance matrix C 1 as: 



— */» 



(3.8) 



where C" 1 e ^ x ^ l \ 

—in-\) 



where C 1 e^ _1) and where C 1 eR. 
— *j* ~ M 
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The variable u stands for 



v T C X 

■V-l) 



The conditional density j > (y„ly (rt _ 1) ) is therefore the prob- 
ability density of a Gaussian normal distribution with a mean 
of u and variance of 10 



l 



It must hold that C^O in order for the above variance to 
exist. This is ensured on the basis of the following argumen- 
tation. Let C and C" 1 have the same characteristic directions 
and inverse~eigenvalues. An eigenvalue 0 of the matrix C" 1 20 
would therefore have an infinite variance of the random vari- 
able vector (AB FBA X1), . . . ,&B FB Jn)) as a consequence. It 
may therefore be presupposed that all the eigenvalues of C 1 
are nonvanishing. Since the eigenvalues of C -1 are non-nega- 
tive, in any case, it therefore holds that: The matrix C" 1 is 25 
symmetrical and positively defined. It is possible by renam- 
ing the coordinate axes to bring this matrix from the form 
(3.8) to the following form: 



T 1 :=: 



r 1 



This matrix is likewise symmetrical and positively defined 
per constructionem. In accordance with the Sylvester crite- 
rion for symmetrical and positively defined matrices, it fol- 
lows therefrom that (C _1 ) n =C,7^>0, and the assertion is 
demonstrated. 

The method according to the invention reduces a simula- 
tion of 



noise sources to the generation of Gaussian-distributed ran- 
dom numbers. 

The inverted covariance matrix C~~ l (an nxn matrix) is 
required in order to generate a random number y rt that is 
correlated in the way required with an already generated 
sequence y (n _ 0 . Strictly speaking, all that is required is 
knowledge of the nth row of this matrix, that is to say knowl- 
edge of 



((rr 1 ) r ,£T 1 ). 



As is to be read from formula (3.6), the covariance matrix C 
depends on the decomposition of the simulation interval [0,0 
into disjoint subintervals (step sizes) [t,_, , tj. In particular, the 
last column of C depends on t„ (doing so identically with the 65 
last row because of the symmetry of C), and thus on the 
current step size At„=t„— t„_ T . 



The left-hand upper (n-l)x(n-l) submatrix U of the nxn 
covariance matrix C is precisely the covariance matrix for a 
random number sequence of length n-1. This covariance 

(3.9) matrix had already to be determined and inverted for calcu- 
lating y(n-l) or for calculating the last element (y^). hi 
order to accelerate the method, it is therefore possible to have 
recourse to an incremental method for matrix inversion, for 
example by use of the Schur complement. 

The invention is also implemented in a method for simu- 
lating a technical system that is subject to 1/f noise. In this 
case, random numbers that have been determined with the aid 
of a method according to the invention are used in modeling 
and/or fixing the variables present on the input channels of the 
system. 

15 Also provided is a computer system and/or a computer 
program for determining sequences of random numbers of a 
1/f noise, or for executing the other methods according to the 
invention. The invention is also implemented in a data 
medium having such a computer program. Furthermore, the 
invention is implemented in a method in which a computer 
program according to the invention is downloaded from an 
electronic data network, such as from the Internet, for 
example, on a computer connected to the data network. 

Other features which are considered as characteristic for 
the invention are set forth in the appended claims. 

Although the invention is illustrated and described herein 
as embodied in a method for on-demand generation of indi- 
vidual random numbers of a sequence of random numbers of 
1/f noise, it is nevertheless not intended to be limited to the 
30 details shown, since various modifications and structural 

(3.10) changes may be made therein without departing from the 
spirit of the invention and within the scope and range of 
equivalents of the claims. 

The construction and method of operation of the invention, 
however, together with additional objects and advantages 
thereof will be best understood from the following descrip- 
tion of specific embodiments when read in connection with 
the accompanying drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a block diagram of a technical system to be 
simulated; 

FIG. 2 is a structogram for determining sequences of ran- 
dom numbers of 1/f noise; 

FIGS. 3A-3F show calculation examples for a first simu- 
lation time step; 

FIGS. 4A-4F show calculation examples for a second 
simulation time step; and 

FIGS. 5A-5F show calculation examples for a third simu- 
lation time step. 

DESCRIPTION OF THE PREFERRED 
EMBODIMENTS 

Referring now to the figures of the drawing in detail and 
first, particularly, to FIG. 1 thereof, there is shown a sche- 
matic diagram of a noisy system that is to be simulated. 
The system is described by a system model 1 that is 
60 indexed as a box and describes the system behavior. The 
system behavior results from input channels 2, which are also 
designated as vector INPUT and from output channels 3, 
which are also designated as OUTPUT. Furthermore, a sys- 
tem-dictated noise is provided, which is present on noise 
input channels 4 and which is also designated as a vector or as 
a matrix NOISE. The matrix NOISE is present when the noise 
is taken into account with a plurality of channels, each col- 
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limn of the matrix NOISE containing a vector of noise values n from 1/f-distributed random numbers. The vector y* serves 
which are present on a noise input channel. for simulating one of the noise sources for the first n simula- 
te noise on the noise input channels 4 is preferably inter- tion time steps, 
preted as a noise-dictated alteration of the system model 1. Eac h dement y k of the nth row of NOISE is then deter- 

The behavior of the input channels 2 and of the output 5 on me basis of equa tions (3.7M3.9), from the last 

channels 3 can be described by a system of differential equa- random number x^ of the vector x and the variables jx x and o, 

tions or by a system of differential-algebraic equations, so tQ be recise according to the following rule: 
that reliable predictions of the system behavior are possible. 

For each time step of the simulation of the system shown in y^*^. 

FIG. 1, a vector OUTPUT of the output channels 3 is calcu- 10 embodiments reproduc- 

lated for a vector INPUT present on the input channels 2 and . nua.jAw^ cp f J 

for a vector NOISE present on the noise input channels 4. *g concrete results of calculation^ 

For thesimulation over arelatively long period of time, the The value of the spectral value p is always adopted in this 

vectors INPUT OUTPUT, NOISE are expediently specified case as 0.5. The value of the intensity const is arbitrarily 

as a matrix, a respective column k of the relevant matrix is adopted as 1.0. Three random numbers are simultaneously 

containing the values of the corresponding time series of the processed in each case in accordance with the simulation of 

relevant INPUT, OUTPUT, NOISE. three noise sources simultaneously acting, on separate chan- 

FIG. 2 illustrates how a respective vector y* is attained nels, on the system to be simulated and which are disposed in 

which forms a column k of the matrix NOISE for the noise a vector y^ k being an integral value from 1 to 3. 

input channels 4 of the system model 1 . Each vector y* serves 20 YIGS . 3 A-3F show a calculation example for a first simu- 

for simulating a noise source. lation time step [to,ti]=[0,0.5]. 

In a first step, a desired spectral value % and an intensity piG 3A shows ^ C0V ariance matrix C of dimension lxl 

constant const are fixed. Furthermore, a counter n of the f or generating a 1/f-distributed random number in the case of 

current simulation time interval is set to 0. The following ^ e simu i at i on step s ize. q represents here only a scalar with 

sequence of computation steps is then carried out succes- 25 ^ Q 7Qj be cmsQ C(l,l) — that is to say with i=j=l— 

sively for each simulation time step. results upon applying equation (3 .6) as 

First, the current simulation time step is fixed. Equivalently o ^ 

to this, it is also possible to fix the end of the current Simula- 1 -°<- |r ^;^-^ 

tion time step, thereby producing the next point in time under <mI " ya °* Q - 5 

consideration. - 30 fig. 3B shows the inverse of the covariance matrix C from 

Afterward, the counter n of the current simulation time step FIG 3Aj some thing which was performed here by use of a 

is incremented by one. Cholesky decomposition not illustrated in more detail here. A 

The covariance matrix C of dimension (nxn) is subse- check 0 f (C C~ 1 H0.707106 . . . 0.707106 . . . _1 ) yields the 

quently determined according to equation (3.6). ^ correct value 1, and this illustrates the correctness of the value 

This is followed by the step of inverting the matrix C, for for C(l,l). 

example by use of a Cholesky decomposition. In order to raise FIG 3C ghows a variable o for the first simulation step n=l . 

the efficiency, it is also possible in this case to appeal to the ^ i s y^ed fr om the equation 

inverse matrix of the previous step, for example by using 

Schur complement techniques. 40 0-^(1/0. 7071 06 . . . ), 

Next, the variable a is calculated from the formula " . .' /in , , nfl t!l _ ~ lfnr , pn t 

sqrt denoting the square root and e ( 1 , 1 ) denoting the element 
<5=sqrt(\/e(n,n)) 0.707106 . . . , indexed by (1,1), of the inverted covariance 

matrix (3 — ^ 

where sqrt denotes the "square root" function and where - ' *w«i1iipq v * x of a CO lV-normallv 

C" wbchismdexedby^. plated in each case, devalues form the first element of 

hiaddition, a value of a (0 1 ^-normally distritoedrandom distributed random num- 

vanable X k is extracted and the vector x* of to normally ^ ^ numbers have the expectation of 

distributed random numbers is thereby supplemented. The ±uc 
extracted random number has the expectation of 0 and the 50 - u ana me vanance 01 ■ 

variance of 1 . This step is carriedout foreachnoise source to FIG. 3E shows three vanables for each of the Jree noise 
h P simulated sources to be simulated. The variable is yielded m accor- 

Furthermore, a variable u fc is formed. It is fonned from the dance with formula (3 .9) from ^^mM of 
first (n-1) components of toe nth row of the inverted covari- the nth row of the inverted covariance matrix C as weU as 
aTceltrix C- and from the sequence of (n-1) 1/f-distrib- 55 from the sequence of (n-1) l^^^^S 
uted random^numbers which were calculated for the preced- that were calculated for me preceding (n-1) s^ulauon time 
\Z (n-\ ^ simulation time steos steps. In the first simulation step, these two vectors each have 

^SSSoiisoiaccordancewtofor. the lengthO. Tt therefore foDows for all the vanables u fc in the 
mula (3.9). This step is carried out for each noise source k to first simulation step that: u*=0. 

be simulated. «o FIG.SFshowsthreevectorsyAoflengthlofl/f-distnbuted 

Finally that element of the matrix NOISE is calculated random numbers that simulate the behavior of three 1/f-dis- 
whose column index k specifies the noise source to be simu- tributed noise sources of the first simulation time : step LO, 
lated and whose row index is equal to n. The current Simula- t,]=[0, 0.5]. The matrix NOISE is yielded from the three 
tion time step is denoted by this. The currently calculated vectorsy^.Thevaluekisanmtegralvalueoffromn 
element r(k,n) of the matrix NOISE represents a random 65 case.Eachelementy.ofmeflrstrowofNOBEisdetermme^ 
number which, together with the superjacent (n-1) elements on the basis of equations (3.7M3.9) according to the lollow- 
of the same column k of NOISE, forms a vector y* of length ing rule from the last random number x* of Ihe associated 
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vector x and the variables u^d a,. For example, the first 
element y^n^ ) of the first vector y, is calculated as follows: 



Vi(?i=l)=x 1 («=l)* o+m=-0.35 . . . 
*0.84 . . . 40.00 . . . =-0.30 . . . 

y 2 (n=l) and ¥3(n=l) are calculated by analogy herewith. 

With the aid o f FIGS . 4A to 4F, there is shown a calculation 
example of a second simulation time step [tj, t 2 ]=[0.5, 0.75]. 
The value n for the second simulation time step is always 
equal to 2. 

FIG. 4A shows the covariance matrix C of dimension 
(nxn)=2x2, which is required for generating one further ran- 
dom number each per noise source. Together with the result in 
accordance with FIG. 3F, the random number thus newly 
generated forms a vector y* of length 2 from 1/f-distributed 
random numbers. A vector y^ is generated in this case for each 
noise source. The covariance matrix C is determined in this 
case according to equation (3.6). 

For example, this is carried out at the element C(2, 1 ) — that 
is to say with i=2 and j=l. From applying equation (3.6), 
C(2,l) is yielded as 



12 



FIG. 4E shows three variables \i k that were calculated from 
the inverted covariance matrix C" 1 in accordance with step 4B 
and from the three random numbers in accordance with step 
3F. For each noise source to be simulated, the variable \i k is 

5 calculated from the (n-1) first components of.the nth row of 
the inverted covariance matrix C" 1 and from the sequence of 
a quantity (n- 1 ) of 1 /f-distributed random numbers that were 
calculated in accordance with formula (3.9) for the preceding 

10 (n- 1 ) simulation time steps . In the second simulation step, the 
variable u fc is thus calculated from the first component of the 
second row of CT 1 as well as from the first component of the 
vector Tnis is carried out, for example, with the aid of the 

l5 value m: 



20 



-0.30...— 1.15.. 



4.79.. 



= -0.07.. 



25 



iri-/2-il°- 5+1 -ki-i-'2-il a5+I ) = 
H0.5-0.75|°- 5+1 + |0-0.75|°- 5+1 + 
|0.5-0.5|° 5+1 -|0-0.5| a5+l )= 30 
-0.125 + 0.6495... + 0-0.3535... = 0.1709... 



35 



FIG. 4B shows the inverse of the covariance matrix C from 
FIG. 4A. A check (not illustrated here) of the condition (C 
C" 1 ) yields a matrix of dimension 2x2, in the case of which 
the elements indexed by (1,1) and (2,2) are equal to 1. The 
other elements have the value 0. 

FIG. 4C shows a variable a that is calculated from the 40 
inverted covariance matrix CT 1 of step 4B. The variable o is 
yielded as 

=sqrt(V4. 



<j=sqrt(\/e(n,n))nsqrt(\/e(2 y 2))= 
79 . . . )=0.45 . . . ; 



45 



50 



sqrt denoting the square root, ande(2,2) denoting the element, 
indexed by (2,2) of the inverted covariance matrix C" 1 from 
FIG. 4B. 

FIG. 4D shows three vectors x* of independent (0,l)-nor- 
mally distributed random numbers, the vectors x* each having 
a length of 2. A (0, 1 )-normally distributed random variable \ k 
is extracted per noise source to be simulated. The extracted 
random number in each case has the expectation of 0 and the 
variance of 1. The vectors x* of the normally distributed 55 
random numbers from FIG. 3D are thereby supplemented 
such that the vectors x* of the normally distributed random 
numbers from FIG. 4D result. 



FIG. 4F shows three vectors y* of length 2 with 1/f-distrib- 
uted random numbers that simulate the behavior of three 
1/f-distributed noise sources for the second simulation time 
step ft , t 2 ]=[0.5, 0.75]. The matrixNOISE is yielded from the 
three vectors y*. The value k is an integral value of from 1 to 
3 in this case. Each element y k of the second row of NOISE is 
determined on the basis of equations (37H39) according to 
the following rule from the last random number x k of the 
associated vector x and the variables u* and o. By way of 
example, the second element y^n=2) of the first vector ^ is 
calculated below: 



ifi("=2)=x 1 (»=2)*a+^ l = 
-0.07 .. .==0.10. . 



0.39 ... * 0.45 . . . 



With the aid of FIGS. 5A-5F, there is shown a calculation 
example for a third simulation time step [t 2 , t 3 ]=[0.75, 1 .25]. 
The value n is always equal to 3 during the third simulation 
time step. 

FIG. 5A shows the covariance matrix C of dimension 
(nxn)=3x3, which is required for generating one further ran- 
dom number each per noise source. Together with the result in 
accordance with FIG. 4F, the random- number thus- newly 
generated forms a vector y^ of length 3 from 1/f-distributed 
random numbers. A vectory^is generated in this case for each 
noise source. The covariance matrix C is determined in this 
case according to equation (3.6). 

For example, this is carried out at the element C(3 , 1 )— that 
is to say with i=3 and j=l. From applying equation (3.6), 
C(3,l)is yielded as 



1.0 (-1/1 +1*1-1 -f 3 l°- 5+1 +1/1 -'3-ll°- 5+1 "I '1-1 -'3-1 l°- 5+1 ) ■ 

= (-|0.5-1.25[°- 5+1 +|0-1.25| 05+1 + jO.5-0.751^ 1 -|0-0.75(° 5+1 ) 
= -0.6495... +1.3975... +0.125 - 0.6495... =0.22... . 
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FIG. 5B shows the inverse C" 1 of the covariance matrix C 
from FIG. 5A. A check (not illustrated here) of the condition 
(C C~ l ) yields a matrix of dimension 3x3, in the case of which 
theelements indexed by (1,1), (2,2) and (3,3) are equal to 1 . 
The other elements have the value 0. 

FIG. 5C shows a variable a that is calculated from the 
inverted covariance matrix C~ l of step 5B. The variable a is 
yielded as 

<s=sqrt{\ /e(n,n)y=sqrt(\/eQ ,3)y=sqrt 
(1/1.75 . . . >=0.75 . . . ; 

sqrt denoting the square root, and e(3 ,3) denoting the element, 
indexed by (3,3) of the inverted covariance matrix C" 1 from 
FIG. 5B. 

FIG. 5D shows three vectors x* of independent (0,1) nor- 
mally distributed random numbers, the vectors x* each having 
a length of 3 . A (0, 1 ) normally distributed random variable x k 
is extracted per noise source to be simulated. The extracted 
random number in each case has the expectation of 0 and the 20 
variance of 1. The vectors x* of the normally distributed 
random numbers from FIG. 4D are thereby supplemented 
such that the vectors x* of the normally distributed random 
numbers from FIG. 5D result. 

FIG. 5E shows three variables \i k that were calculated from 25 
the inverted covariance matrix CT 1 in accordance with step 5b 
and from the three random numbers in accordance with step 
4/ For each noise source to be simulated, the variable u fc is 
calculated from the (n-1) first components of the nth row of 
the inverted covariance matrix C" 1 and from the sequence of 30 
(n- 1 ) 1 /f-distributed random numbers that were calculated in 
accordance with formula (3.9) for the preceding (n-1) simu- 
lation time steps. In the second simulation step, the variable \i k 
is thus calculated from the first two components of the third 
row of C" 1 as well as from the first two components of the 35 
vector y^. This is carried out, for example, with the aid of the 
value \i Y : 



Ml 



= n.n 

-0.30... -0.31... + 0.10... -0.98 



1.75... 



= -0.00.... 



40 
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FIG. 5F shows three vectors y* of length 3 with 1/f-distrib- 
uted random numbers that simulate the behavior of three 
1 /f-distributed noise sources for the third simulation time step 50 
[t 2 , t 3 ]=[0.75, 1.25]. The matrix NOISE is yielded from the 
three vectors y*. The value k is an integral value of from 1 to 
3 in this case. Each element y k of the third row of NOISE is 
determined on the basis of equations (3.7H3.9) according to 
the following rule from the last random number x*(n=3) of the 55 
associated vector x and the variables \i k and o. By way of 
example, the third element y_,(n=3) of the first vector y, is 
calculated below: 



j/ l (n=3)=jc 1 («=3) *a-ni 1 < 
+0.00 . . . =0.67 . 



=-0.90 ... * 0.75 . 



60 



for a third simulation interval. In this case, after accurate 
numerical calculation, all values are truncated after the sec- 
ond decimal place following the point in order to be able to 
reproduce them more effectively. Consequently, when com- 
putationally reaccomplishing the exemplary embodiments, it 
is necessary to calculate further not with the intermediate 
values shown in the figures, but with the exact intermediate 
values, in order to reach the specified y vectors starting from 
the specified x vectors. 

Vectors of (0, 1) normaUy distributed random variables are 
shown in FIGS. 3C, 4C and 5C. In this case, a random variable 
represents a noise source in each case. For the sake of sim- 
plicity, it is not illustrated here how such random numbers are 
reached with the aid of the expectation of 0 and the variance 
of 1 . This is familiar to the person skilled in the art. 

We claim:. 

1. A computer implemented method for numeric simula- 
tion of an electric circuit influenced by noise, comprising the 
steps of: 

numerically simulating the electric circuit using a model 
including input channels, noise input channels and out- 
put channels, the behavior of the input channels and of 
the output channels being described by a system of dif- 
ferential equations or a by a system of differential alge- 
braic equations; 
in the numerical simulation step, calculating an output 
vector for an input vector present on the input channels 
and a noise vector y of 1/f distributed random numbers 
present on the noise input channel, wherein generation 
of the noise vector y includes the steps of: 
determining a desired spectral value p; 
determining a number of the random numbers of the 1/f 

noise to be generated; 
determining an intensity constant const; 
setting a starting value for a running variable n; 
performing a loop-type repetition until a desired number of 
elements y(n) of a vector y of length n is calculated from 
1 /f-distributed random numbers, by the steps of: 
increasing a current value of the running variable n by 1 ; 
setting a simulation time step [t^ ; t J ; 
determining elements C, y of a covariance matrix C of 
dimension (nxn) according to: 

17=1 n 

determining an inverted covariance matrix C" 1 by inverting 

the covariance matrix C ; 
detennining a variable o in accordance with 

0=sqrt{\/e{n,n) ), 

where sqrt denotes a square root function, and e(n,n) 
denotes an element of the inverted covariance matrix 
CT 1 indexed by (n,n); 
determining a (0, l)-normally distributed random number 
which forms an nth component of a vector x of length n; 
forming a variable u from first (n-1) components of an nth 
row of the inverted covariance matrix C" 1 and (n-1) 
elements o f the vector y calculated for a preceding (n- 1 ) 
simulation time step, according to: 



The following conditions remain to be observed for the 
purpose of concrete execution of the calculation examples 
shown. 

The numerical values shown in FIGS. 3, 4 and 5 reproduce 65 
intermediate and final results of the computing steps, 
described with reference to FIG . 2, for a first, for a second and 
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where y (n _D denotes first (n-1) elements of the vector y, 
CrJ denotes the first (n-1) components of the nth row of 
the inverted co variance matrix C~\ and denotes a 
component of the inverted covariance matrix C~ 
indexed by (n,n); and 5 

calculating an element y (n) of the vector y of length n from 
the 1/f-distributed random numbers, according to: 

outputting the at least one sequence of random numbers of 1/f l o 
noise. 

2. A method for simulating a technical system subject to 1/f 
noise, which comprises the steps of: 

determining random numbers according to claim 1; and 
using the random numbers for modeling variables present 
on input channels of the technical system. . 

3. A computer system, comprising: 

processor programmed for executing the method for deter- 
mining the sequences of random numbers of the 1/f 
noise according to claim 2. 

4. A computer ninning a computer program, comprising: 
computer-executable instructions for carrying out the 

method according to claim 1 for detemiining the 
sequences of random numbers of the 1/f noise. 

5. A computer running the instructions stored on a com- 
puter-readable data medium having the computer-executable 
instructions according to claim 4. 

6. A downloading method, which comprises the step of: 
downloading the computer program according to claim 4 

from an electronic data network onto a computer con- 
nected to the electronic data network. 

7. The method according to claim 6, which further com- 
prises using the Internet as the electronic data network. 

8. A computer system, comprising: 
processor programmed for executing the method for deter- 

mining the sequences of random numbers of the 1/f 
noise according to claim 1. 

9. Amethod for simulating a technical system subject to 1/f 
noise, which comprises the steps of: 

determining random numbers according to claim 1; and 
using the random numbers for fixing variables present on 

input channels of the technical system; and 
simulating the technical system and outputting the result of 

the simulation. 

10. A computer implemented method for numeric simula- 
tion of an electric circuit influenced by noise, comprising the 
steps of: 

numerically simulating the electric circuit using a model 
*~" including input channels, noise input channels and out- 
put channels, the behavior of the input channels and of 
the output channels being described by a system of dif- 
ferential equations or by a system of differential alge- 
braic equations; 

in the numerical simulation step, calculating an output 
vector for an input vector present on the input channels 
and a noise vector y of 1/f distributed random numbers 
present on the noise input channel, wherein generation 
of the noise vector y includes the steps of: 

determining a desired spectral value P; 

determining a number of the random numbers of the 1/f 
noise to be generated; 

determining an intensity constant const; 

setting a starting value for a running variable n; 

calculating q sequences of the random numbers of the 1/f 
noise simultaneously, by performing loop-type repea- 
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tions until a desired number of elements y kj „ of a vector 
y of length n is calculated from 1/f-distributed random 
numbers, by the steps of: 

increasing a current value of the running variable n by 1 ; 
setting a simulation time step [t„_ x ; t J ; 
determining elements of a covariance matrix C of 
dimension (nxn) according to: 

determining an inverted covariance matrix C" 1 by inverting 

the covariance matrix C ; 
detennining a variable a in accordance with 

G=sqrt (\/e(rt,n) ), 

where sqrt denotes a square root function, and e(n,n) 
denotes an element of the inverted covariance matrix 
C" 1 indexed by (n,n); 

determining a quantity q of (0, 1 )-normally distributed ran- 
dom numbers x k>n which form a respective last compo- 
nent of vectors x^ of length n, where k=l , . . . , q, 

forming q variables according to: 



c- 1 



where y (n ^ x * denotes first (n-1) elements of the vector y* 

that were' calculated for a preceding simulation time 

step, Cr^ 1 denotes the first (n-1) 
components of the nth row of the inverted covariance 

matrix C" 1 , and C^ n denotes the element of the inverted 

covariance matrix C" 1 indexed by (n,n), where 

k=l, . . . , q; and 
calculating q elements y kt „ which form a respective nth 

component of the vector y* of length n from 1 /f-distrib- 

uted random numbers, according to: 

where k=l, . . . , q; 

outputting at least one of the q sequences of random num- 
bers of 1/f noise. 

11. A method for simulating a technical system subject to 
1/f noise, which comprises the steps of: 

determining random numbers according to claim 10; and 

using the random numbers for modeling variables present 
on input channels of the technical system. 

12 A computer system, comprising : 

processor programmed for executing the method for deter- 
mining the sequences of random numbers of the 1/f 
noise according to claim 11. 

13. A computer running a computer program, comprising: 
computer-executable instructions for carrying out the 

method according to claim 10 for determining the 
sequences of random numbers of the 1/f noise. 

14. A computer r unnin g the instructions stored on a com- 
puter-readable data medium having the computer-executable 
instructions according to claim 13. 

15. A downloading method, which comprises the step of: 
downloading the computer program according to claim 13 

from an electronic data network onto a computer con- 
nected to the electronic data network. 

16. The method according to claim 15, which further com- 
prises using the Internet as the electronic data network. 
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17. A computer system, comprising: 

processor programmed for executing the method for deter- 
mining the sequences of random numbers of the 1/f 
noise according to claim 10. 

18. A method for simulating a technical system subject to 
1/f noise, which comprises the steps of: 
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determining random numbers according to claim 10; and 
using the random numbers for fixing variables present on 

input channels of the technical system; and 
simulating the technical system and outputting the result of 

the simulation. 



